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Highlight 
e Two pheasant coronaviruses (PhCoVs) were isolated in 2017 in China. 
e The two PhCoVs were genetically similar to IBV. 
e Pathogenicity, replication, and shedding of PhCoV were obvious different 
when infected chickens and pheasants. 
e PhCoVs isolated from different outbreaks may have evolved independently 


from IBVs by adaption in pheasants. 


ABSTRACT 

Two viruses were isolated in 2017 from commercial pheasants with severe 
clinical signs and mortality in Shandong and Anhui provinces, China, respectively. We 
examined the pathogenic effects of the viruses in chicken embryos and the size and 
morphology of the virus particles, performed phylogenetic analysis based on the S1 
gene and complete genomic sequences, and examined the antibody responses against 
infectious bronchitis virus (IBV). The results suggested that the viruses [0623/17 and 
[0710/17 were avian coronaviruses and were identified as pheasant coronaviruses 
(PhCoV), with greatest similarity to IBV. Further investigations of the antigenicity, 
complete genome organization, substitutions in multiple genes, and viral 
pathogenicity, replication, and shedding in chickens and pheasants showed obvious 
differences between PhCoV and IBV in terms of antigenicity, and viral pathogenicity, 
replication, and shedding in chickens and pheasants. The close genetic relationship, 
but obvious differences between PhCoVs and IBVs suggested the IBVs could be the 
ancestors of PhCoVs, and that PhCoVs isolated from different outbreaks may have 
evolved independently from IBVs circulating in the specific region by adaption in 
pheasants. This hypothesis was supported by analysis of the S1 gene fragments of the 
two PhCoVs isolated in the current study, as well as PhCoVs isolated in the UK and 
selected IBV strains. Such analyses indicated different evolution patterns and different 
tissue tropisms between PhCoVs isolated in different outbreaks. Further studies are 


needed to confirm this hypothesis by studying the complete genomic sequences of 


PhCoVs from different outbreaks and the pathogenicity of IBVs in pheasants to 


compare and clarify the relationships between PhCoVs and IBVs. 
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1. Introduction 

Coronaviruses (CoVs) are enveloped, positive-sense RNA viruses found in a wide 
range of animals, including mammals and birds. CoVs exhibit marked tropism for 
epithelial cells in the respiratory, digestive, and urogenital tracts, and are responsible 
for a diverse spectrum of enteric, hepatic, neurological, and respiratory illnesses of 
differing severities. CoVs, classified in the subfamily Coronavirinae, family 
Coronaviridae, order Nidovirales, are now recognized as emerging viruses with a 
propensity to cross into new host species (Leopardi et al., 2018). CoVs are subdivided 
into four genera on the basis of genotypic and serological characterization: alpha and 
beta coronaviruses, which are found mainly in mammals; gamma coronaviruses, 
detected in birds and marine mammals; and delta coronaviruses found mainly in birds 
and swine (de Groot, 2012). The 5' two-thirds of the CoV genome encodes proteins 
involved in viral RNA synthesis. Most of these proteins are characteristically encoded 
by the partially-overlapping open reading frames (ORFs), ORFla and ORF 1b, and are 
translated as the polyproteins ppla and pplab, which are then processed by virus- 


encoded proteinases into 15 or 16 non-structural proteins (nsps) (Ziebuhr, 2005). The 


remaining 3' one-third of the genome encodes virus structural proteins, including the 
spike (S) glycoprotein, membrane (M) glycoprotein, small envelope (E) protein, and 
phosphorylated nucleocapsid (N) proteins. The CoV S protein can be divided into an 
amino-proximal half (S1 subunit or domain) containing the receptor-binding domain, 
and a carboxyl-proximal half (S2 subunit or domain) containing elements involved in 
membrane fusion (Lewis et al., 2015). The S protein is an important target for T cell 
responses and is the major inducer of virus-neutralizing antibodies, which are elicited 
by epitopes located mostly in the molecule’s S1 domain (Reguera et al., 2012; Satoh 
et al., 2011). The S1 and S2 domains in some CoVs, e.g. avian CoVs (AvCoYVs), are 
cleaved by a cellular furin-like enzyme (de Haan et al., 2014). 

Most viruses in the genus Gammacoronavirus are isolated from birds, including 
Galliformes (chicken, turkey, quail, guinea fowl, pheasant, and peafowl) (Brown et 
al., 2016; Cavanagh et al., 2002; Guy, 2008; Liu et al., 2005; Sun et al., 2007; Torres 
et al., 2017), Anseriformes (duck, goose, teal, swan, and pintail) (Liu et al., 2005; 
Papineau et al., 2019), Columbiformes (pigeon) (Jonassen et al., 2005), 
Pelecaniformes (spoonbill and heron), Suliformes (cormorant), Charadriiformes (red 
knot, oystercatcher, and black-headed gull), and Passeriformes (bulbul) (Chu et al., 
2011), although they have also been identified in the beluga whale 
(Mihindukulasuriya et al., 2008) and in bottlenose dolphins (Woo et al., 2014). 
Infectious bronchitis virus (IBV) and turkey coronavirus (TCoV) are among the most 
economically important AvVCoVs. IBV was the first recognized CoV in the United 


States in 1931 (Schalk et al., 1931). IBV is traditionally considered to be a host- 


specific pathogen in chickens, responsible for respiratory, renal, and genital diseases, 
with heavy economic consequences worldwide (Cavanagh, 2007). In contrast to alpha 
and beta coronaviruses, which only occur as one or two different serotypes, [BV has 
many different serotypes, genotypes, lineages, and variants (Valastro et al., 2016). 
Furthermore, new IBV genotypes/lineages and variants are continuing to emerge due 
to its high rate of evolution. This is expressed as an accelerated rate at which viable 
point mutations, nucleotide insertions, or deletions accumulate in the genome, 
especially in the S1 gene, associated with a lack of proofreading viral polymerase in 
IBVs (Cavanagh, 2007), as well as the high rate of virus replication. In addition, 
evidence suggests that some IBV strains may have arisen by genetic recombination 
(Cavanagh et al., 1992; Kusters et al., 1989; 1990), or by a combination of mutation 
and recombination (Jia et al., 1995; Chen et al., 2016; Jiang et al., 2017). TCoV, 
initially found in the 1970s, is associated with the enteric disease known as 
transmissible enteritis, coronaviral enteritis of turkeys, or bluecomb (Guy, 2008). 
TCoV infection was also recently found to be associated with a syndrome comprising 
several intestinal disorders, usually occurring in turkeys within the first 3 weeks of 
life and referred to as poult enteritis complex (Barnes et al., 2000). TCoV emergence 
was proposed to be the result of recombination events involving IBVs and an 
unidentified CoV, which donated an S gene encoding a protein of low amino acid 
identity to those of IBVs (Hughes, 2011; Jackwood et al., 2010). These 
recombinations resulted in a host shift from chickens to turkeys, and in a change in 


tissue tropism of the virus from the upper respiratory tract to the intestine. 


Outbreaks initially associated with CoV infection in pheasants were first reported 
in the UK in 1980 (Spackman and Cameron, 1983). Affected birds exhibited reduced 
egg production among breeding hens, accompanied by loss of shell pigmentation and 
some abnormal shell quality, with later outbreaks of respiratory signs among affected 
pheasants. However, there was no evidence of kidney damage in the affected birds. 
Antibodies against IBV were detected by the hemagglutination inhibition (HI) test in 
80% of tested pheasants and a CoV was isolated that proved to be identical to the 
agent responsible for infectious bronchitis (IB) in chickens. Two later outbreaks of 
nephritis occurred in the UK in 1983 and 1984, resulting in the deaths of 450 out of 
1000 8-week-old pheasants and 15 out of 160 adult pheasants, respectively, from 
which CoVs were isolated using 8—9-day-old embryonating eggs (Lister et al., 1985). 
However, intranasal instillation of specific-pathogen-free (SPF) chickens with the 
virus failed to produce any clinical signs of disease. A similar outbreak of nephritis in 
1994 resulting in the deaths of over 1000 breeding pheasants out of a total stock of 
approximately 7000 birds in the UK was considered to be linked to avian IBV (Gough 
et al., 1996). Sneezing and reduced egg production and hatchability were also 
reported, but the egg quality was unaffected. Pheasant coronavirus (PhCoV)- 
associated nephritis and respiratory disease were also subsequently reported in the UK 
(Pennycott, 2000; Cavanagh et al., 2002). Partial genomic sequence analysis showed 
that the gene sequences of the pheasant viruses differed from those of IBV to a similar 
extent as the sequence of one serotype of IBV differs from another (Cavanagh et al., 


2002). The above results showed that the CoVs involved in these pheasant infections 


were closely related to avian IBV; however, the clinical signs and mortality associated 
with CoV infections in pheasants and the isolated viruses differed from each other and 
from IBVs. Further studies using molecular techniques are therefore required to 
determine the complete genomic sequences of CoVs from different outbreaks in 
pheasants, to compare and fully elucidate their origins and antigenic relationships 
with strains of IBV. 

In this study, we isolated two CoVs from pheasants with nephritis and high 
mortality for the first time in China. We aimed to determine the full genome 
sequences of the two isolates and compare them with IBVs, evaluate the 
pathogenicity of the isolates in pheasants and chickens, and investigate the cross- 
antigenicity with IBVs, to clarify the origin of PhCoVs and their relationships with 


strains of IBV. 


2. Materials and methods 
2.1. Chicken embryos, chickens, and pheasants 

White leghorn chickens were hatched from SPF embryonated eggs (Harbin 
Veterinary Research Institute, Chinese Academy of Agricultural Sciences), maintained 
in Hosfall-type isolators, and provided with food and water ad libitum. Chicken 
embryonated eggs were used for primary isolation of viruses from the diseased 
pheasants, viral seed stock preparation and titration, virus-neutralization (VN) tests, 
re-isolation attempts from oropharyngeal and cloacal swabs, and titration of viruses 


from selected tissues from chickens and pheasants in the challenge trials. SPF 


chickens were used for antisera production and challenge trials. One-day-old 
commercial pheasants (Phasianus colchicus) were purchased from a pheasant- 
producing company in Shandong province, China, and used in the challenge trials. All 
experimental procedures and animal care were approved by the Ethical and Animal 
Welfare Committee of Heilongjiang province, China (License nos. HS Y-[ACUC- 
2017-167, -168, -169, and -170). 
2.2. Clinical information and sample collection 

Several new episodes of serious respiratory diseases occurred in commercial 
pheasant flocks in China between September 2016 and June 2017. The outbreaks 
were characterized by clinical signs compatible with IB in chickens, including 
respiratory signs (all flocks) and kidney lesions (nephritis), leading to economic losses 
in different geographical locations among the major pheasant-farming areas in 
Shandong, Anhui, and Henan provinces in China. The birds showed early signs of 
respiratory disease at around 17 days old, most birds in the flocks showed obvious 
clinical signs at 20-40 days old, and the clinical signs disappeared at around 55 days 
old. Mortality occurred at approximately 20 days old and was generally about 30%, 
but was as high as 50% in some flocks. Gross examination showed severe tracheitis 
and nephritis. Trachea, kidney, and proventriculus samples from two selected farms in 
Shandong and Anhui provinces were sent to our laboratory for virus detection in early 
2017. The two farms contained 13,500 and 7,000 birds, respectively, of which more 
than 4000 and 2500 died, respectively, during the disease outbreaks. 


2.3. Virus detection and isolation 


Tracheal and kidney samples from 10 dead birds from the farm in Shandong 
province, and kidney and proventriculus samples from five dead birds from the farm 
in Anhui province, respectively, were collected. Samples of each farm were pooled 
together, ending up with 2 different pools. The samples tested negative for Newcastle 
disease virus (NDV) (Gohm et al., 2000) and avian influenza virus (AIV) subtypes H5 
and H9 by reverse transcription polymerase chain reaction (RT-PCR) (Chaharaein et 
al., 2009). However, the samples tested positive by RT-PCR using primers initially 
designed to detect the 3' end of the IBV N gene (Liu et al., 2009). 

PhCoVs can replicate and be isolated in embryonated chicken eggs (Cavanagh et 
al., 2002). Therefore, each pooled sample was used to inoculate five 9-day-old 
embryonated eggs (0.2 mL/egg) via the allantoic route. Allantoic fluid was harvested 
from two inoculated eggs for each pooled sample at 72 h post-inoculation, pooled, 
and used for virus detection by RT-PCR and direct electron microscopy (Liu and 
Kong, 2004). The remaining three eggs were incubated at 37°C, candled daily, and 
examined for specific lesions at 7 days post-inoculation (dpi). Further blind passages 
were performed until embryo dwarfing and death were observed at 2—7 dpi. 

Two CoVs were isolated from samples from pheasants from Shandong and Anhui 
provinces, respectively. These PhCoVs were designated yCoV/ph/China/I0623/17 
(10623/17) and yCoV/ph/China/I07 10/17 (10710/17), respectively, according to the 
standard nomenclature adopted for PhCoV isolates in this study (Ducatez et al., 
2016). The two PhCoVs were isolated and propagated by inoculating and passaging in 
the allantoic cavity of 9-day-old SPF embryonated chicken eggs (Liu and Kong, 
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2004). The allantoic fluids from each passage incubated with isolates [0623/17 and 
10710/17 were tested for NDV and AIV by hemagglutination (HA) activity. The 
median titers of isolates 10623/17 and I0710/17 at the final passage (3rd) were 10°° 
and 10°” egg infectious doses (EIDs0), respectively, determined using the method of 
Reed and Muench (1938). 
2.4. Detection of antibody against IBV in commercial pheasants 

Blood samples were collected from six male (A and C) and six layer hen 
pheasants (B and D) from each of two affected commercial pheasant flocks in 
Shandong and Anhui provinces, respectively, approximately 2 months after the 
disappearance of clinical signs. These samples were used for the detection of 
antibodies against IBV using a commercial enzyme-linked immunosorbent assay 
(ELISA; IDEXX Corporation, Westbrook, ME, USA) following the manufacturer’s 
instructions (Liu et al., 2013). Blood samples were also collected from an additional 
six commercial male (E) and six layer pheasants (F) from another flock without 
clinical signs, from the same farms as the affected commercial pheasant flocks in 
Shandong province. None of these birds had received IBV vaccines. 
2.5. RNA extraction, RT-PCR amplification, and sequencing of the complete genome 

The complete genomic RNAs of isolates [0623/17 and I0710/17 were extracted, 
from the respective infective allantoic fluids using a viral RNA extraction kit (Takara 
Bio Inc., Shiga, Japan), according to the manufacturer’s instructions. We initially 
amplified the complete genomes of the two PhCoV isolates using 14 pairs of partly 
overlapping primers, originally used to amplify the complete IBV genome (Liu et al., 
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2013). The sense and antisense primers were also used to amplify the 3’ and 5’ ends of 
the IBV genome (Liu et al., 2013), together with a PrimeScript™ One-Step RT-PCR 
Kit ver. 2 (Takara) and a 3'/5' RACE kit (Takara), according to the manufacturer’s 
instructions. Two pairs of primers and the sense primer for amplifying the 3’ end of 
the virus genome were then redesigned to replace the previous corresponding primers, 
which were unable to amplify the gene fragments of isolates [0623/17 and I0710/17. 
All primers used in the following procedures are available on request from the 
authors. 

The RT-PCR products were resolved by electrophoresis in 1% agarose gels. The 
gel-purified RT-PCR products were subjected to direct sequencing or ligated into the 
PMD-18T cloning vector (Takara), following the manufacturer’s instructions. Three 
to five clones per amplicon were sequenced to determine the consensus sequence for 
any given genomic region. Consensus sequences were assembled from both forward 
and reverse sequences, edited manually, and corrected on the basis of the different 
clones of a fragment. The ORFs and deduced proteins were analyzed using Molecular 
Evolutionary Genetics Analysis 6.0 (Mega 6.0) (Tamura et al., 2013) and 
BioNumerics (Applied Maths, Sint-Martens-Latem, Belgium), by comparison with 
the complete genomic sequences of IBV LX4 (GenBank accession 
numberAY 338732) and H120 (FJ888351) strains, and TCoV TCoV-ATCC 
(EU022526) and TCoV/TX-GL/01 (GQ427174) strains. The sequences reported here 


have been deposited with the public access GenBank database 
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(www.ncbi.nlm.nih.gov/genbank/). The accession numbers of isolates [0623/17 and 
10710/17 are MK423877 and MK423876, respectively. 
2.6. SI gene/fragment analysis 

Phylogenetic trees were constructed by aligning sequences using the Clustal W 
multiple alignment algorithm and the neighbor-joining method with 1000 bootstrap 
replicates, employing the Kimura 2-parameter and nucleotide substitution model. The 
S1 gene sequences of 213 reference AVCOV strains were used for phylogenetic 
analysis, together with those of the two new isolates [0623/17 and I0710/17. The 213 
viruses included 199 reference IBV strains reported by Valastro et al. (2016), three 
IBV GI-28 strains (Chen et al., 2017), three IBV GI-29 strains (Jiang et al., 2018), 
two IBV GVII-1 strains (Ma et al., 2019), two additional IBV GI-19 strains, 
ck/CH/LDL/09 1022 (LDL/091022) (Liu et al., 2013) and ck/CH/LJL/140734 
(LJL/140734) (Zhao et al., 2017) which were used the virus-neutralization (VN) test, 
three TCoVs and one guinea fowl CoV (Brown et al., 2016). 

Nine nucleotide sequences of SI gene fragments of PhCoVs were retrieved from 
the GenBank database. The sequence fragments were then submitted to a BLAST 
search in the NCBI database to identify homologous sequences. Six sequences of IBV 
strains that showed relatively close genetic relationships with the nine PhCoVs were 
selected for phylogenetic analysis and sequence alignment by comparison with our 
two isolates, eight GI-19 strains, and three GI-18 strains. Similarly, 16 nucleotide 
sequences of 3’ UTR fragments of PhCoVs were retrieved from the GenBank 
database and compared phylogenetically with 44 corresponding sequences from 43 
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IBV strains and one duck CoV, selected based on the results of phylogenetic analysis 
of the complete genomes. 
2.7. Comparison of complete genomes 

The complete genomic sequences of PhCoV isolates [0623/17 and [0710/17 were 
compared with those of 84 IBV strains and one duck CoV (as an outgroup) from the 
GenBank database, representing sequences of IBV strains from all continents and 
each lineage. Phylogenetic trees were constructed as described above. In addition, the 
complete genomic sequence of isolate 10623/17 was compared with 39 IBV strains 
and one duck CoV from the GenBank database, which selected on the basis of the 
results of phylogenetic analysis of the complete genome, using Similarity Plot 
(SimPlot) version 3.5.1 (http://sray.med.som.jhmi.edu/SCRopftware/simplot/) to 
analyze possible recombination events. The PhCoV isolate 10710/17 was used as a 
query virus. The window width and step size were set to 1,000 bp and 50 bp, 
respectively. 
2.8. Virus cross-neutralization test 

In addition to the two isolates [0623/17 and [0710/17 in this study, three 
additional IBV strains, LDL/091022, LJL/140734, and H120, were also used in the 
VN test. LDL/091022 and LJL/140734 were selected because they were genetically 
related to isolates 10623/17 and I0710/17 based on the phylogenetic trees using both 
S1 gene/fragment and complete genome sequences. The H120 vaccine was used 
because it was commonly employed in chicken flocks in China, and was also used for 
vaccination of some commercial pheasant flocks in China; the H120 vaccine was not 
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used in the affected pheasant flocks in this study. Relatedness values according to the 
Archetti and Horsfall method (1950) range from O for viruses that are antigenically 
unrelated to 100% for isolates that are identical. This method for assessing the 
antigenic relatedness of IBV means that viruses of the same serotype may have 
relatedness values ranging from 50%—100%, whereas viruses of a different serotype 
typically have relatedness values ranging from 0%—20% (Archetti and Horsfall, 
1950). 

Sera against the five virus strains were prepared as described previously (Guo et 
al., 2014). For VN, sera were serially diluted two-fold with sterile phosphate-buffered 
saline and mixed with 100 EIDso of the virus strains. After incubation for 1 h at 37 °C, 
virus-serum mixtures were inoculated into the allantoic cavity of SPF chicken 
embryos and observed for 7 days. Virus replication was determined by the appearance 
of typical lesions in the inoculated eggs, and by amplification of the IBV N gene from 
allantoic fluid by RT-PCR (Liu et al., 2009). The tests were performed in the presence 
of appropriate controls. The end-point titer of each serum sample and the antigenic 
relatedness values were calculated using the method of Reed and Muench (1938). 

In addition, an aliquot of each serum from male and layer hens from each 
commercial pheasant flock was pooled and tested for PhCoV antibodies by VN, using 
{0623/17 virus. 

2.9. Experimental design 


2.9.1. Chicks 
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A total of 135 1-day-old SPF chicks were divided randomly into nine groups and 
kept in separate isolation units, with 15 chicks per group (Table 1). Chickens in 
groups 1, 4, and 7 were administered 0.1 ml of 10°-° EIDso of the isolate 10623/17 via 
the ocular and nasal routes at 1, 30, and 45 day(s) of age, respectively. Chickens in 
groups 2, 5, and 8 were administered 0.1 ml of 10°-° EIDso of the isolate I0710/17 via 
the ocular and nasal routes at 1, 30, and 45 day(s) of age, respectively. Chickens in 
groups 3, 6, and 9 were administered the same amount of sterilized allantoic fluid at 
the corresponding ages, as negative controls. Five chickens from each inoculated and 
control group were euthanized at 5 dpi, and trachea, lung, proventriculus, kidney, and 
cecal tonsil tissues were separated from the surrounding tissue, removed aseptically, 
and used for virus titration, as described previously (Liu et al., 2009). Trachea and 
kidney tissue samples were used for the detection of IBV antigen using 
immunohistochemistry (IHC) with the monoclonal antibody, 6D10 directed against 
the N protein, as described previously (Han et al., 2013; 2016). The remaining chicks 
were examined daily for signs of infection for another 25 days after inoculation. 
Signs, such as eye irritation and/or scratching of the inoculated eye, lethargy, mild 
coughing, and/or ‘snicking’, were monitored daily for a period of 30 days. Blood 
samples and oropharyngeal and cloacal swabs were collected from the birds in each 
group at 4, 8, 12, 16, 20, 24, and 28 dpi. The blood samples were used for antibody 
detection (ELISA; IDEXX Corporation). The oropharyngeal and cloacal swabs were 
used for virus recovery in 9-day-old SPF chicken eggs (Zhang et al., 2015). 

2.9.2. Commercial pheasants 
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A total of 35 1-day-old commercial pheasants were used in the present study. The 
birds were housed in isolators under negative pressure | week before experimental 
infection (Table 1). Blood samples were collected from all the birds at 5, 10, 15, 20, 
25, 30, 35 and 40 days of age and used for IB V-specific antibody detection (ELISA; 
IDEXX Corporation). The IBV-specific antibodies had not been detected until 40 days 
of age. Five birds were selected randomly and euthanized at 40 days, and examined 
carefully for lesions in each of the tissues and organs. Trachea, lung, kidney, digestive 
tract, liver, oviduct, spleen, and bursa of Fabricius tissue samples were collected and 
examined by histopathology to confirm that the birds used in the study were clinically 
healthy. The samples tested for RT-PCR and PCR were processed differently and 
NDV (Gohm et al., 2000), IBV (Liu et al., 2009), the AIV subtypes H5, H7, and H9 
(Chaharaein et al., 2009) and avian metapneumovirus (Cavanagh et al., 1999) were 
tested by RT-PCR, and fowl adenovirus (Li et al., 2016), Mycoplasma synoviae (MS), 
and Mycoplasma galliscepticum (MG) (Moscoso et al., 2004) were tested by PCR 
using the collected samples. The birds were shown to be free of these pathogen 
infections (data not shown). Blood samples were collected from the remaining 30 
birds. None of them had NDV- or AIV H5-, H7-, or H9-specific HI antibodies at 40 
days of age, and none had antibodies against infectious bursal disease virus, avian 
leukosis virus, reticuloendotheliosis virus, chicken infectious anemia virus, fowl 
adenovirus, MS, or MG according to ELISAs (DEXX Corporation). Then, the flock’s 


health status was observed for another 5 days before the experiments were performed. 
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Thirty commercial pheasants were divided randomly into two groups of 15 birds 
per group. Isolate I0710/17 was selected to infect the pheasants in this study because 
it had the lower rate of virus replication than that of isolate 10623/17 in chickens 
based on the results of infection to SPF chickens. Birds in group 1 were administered 
0.1 ml of 10°-° EIDso of the isolate I0710/17 via the ocular and nasal routes at 45 days 
of age. Birds in group 2 were administered the same amount of sterilized allantoic 
fluid at the same age, as a negative control. Five birds from each group were 
euthanized at 5 dpi, and trachea, lung, proventriculus, kidney, and cecal tonsil tissues 
were collected and used for virus titration. The titers were calculated to include the 
total amount of infectious virus recovered from the entire organ. The tracheas and 
kidneys were also subjected to immunohistochemistry (IHC) for the detection of IBV 

Trachea and kidney samples were used to detect IBV antigen using IHC (Han et 
al., 2013; 2016). The birds were examined daily for signs of infection for another 25 
days after inoculation. Signs were monitored daily for a period of 30 days. Blood 
samples and oropharyngeal and cloacal swabs were collected from the remaining 
birds in each group at 4, 8, 12, 16, 20, 24, and 28 dpi. The blood samples were used 
for antibody detection and the oropharyngeal and cloacal swabs were used for virus 
recovery in 9-day-old SPF chicken eggs (Han et al., 2016), as described above. 

2.10. Statistical analysis 

Data are expressed as the mean + standard deviation. Virus titers were analyzed 
by the Student’s f-test using GraphPad Prism for Windows version 5 (GraphPad 
Software, La Jolla, CA, USA). Differences were considered significant if the p value 
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was <(0.05. 


3. Results 
3.1. IBV-like coronaviruses were isolated from the diseased commercial pheasants 

Tissues from each flock were pooled and inoculated into embryonated SPF 
chicken eggs, as described previously (Liu and Kong, 2004). No obvious lesions were 
observed at 7 dpi in the first and second passages. However, stunting and curling of 
the embryo with urate deposits in the embryonic mesonephros appeared following 
two additional passages in the allantoic cavity of 9-day-old SPF eggs (Gough et al., 
1996). The allantoic fluids from each passage showed no HA activity for chicken red 
cells. Examination of the concentrated allantoic fluids from the third embryo passage 
by electron microscopy revealed CoV-like particles (Supplementary Fig. 1). 

The allantoic fluids tested positive for IBV by RT-PCR using IBV-specific 
primers targeting the 3' end of the N gene, as described previously (Liu et al., 2013). 
To further identify the virus isolates, we sequenced the PCR products of the two 
isolates and used the sequences to search the GenBank database 
(https://blast.ncbi.nim.nih.gov/Blast.cgi) for similar sequences using the BLASTN 
program. The two sequences shared the highest nucleotide identities (90%) with two 
IBV strains, ck/CH/LDL/OSI and ck/CH/LDL/OSII, isolated in China in 2005 and 
belonging to the GI-19 lineage (Zhao et al., 2017). 

All the affected pheasants collected from four commercial flocks 2 months after 
the disappearance of clinical signs were positive for the IB V-specific antibodies by 
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ELISA, in contrast to the negative results for serum from pheasants that had not 
shown clinical signs (Fig. 1). 

Taken together, the embryopathic effect of the virus, the size and morphology of 
the virus particles, the sequence characteristics of the amplified gene fragments, and 
the antibody responses against IBV suggested that the agents 10623/17 and I0710/17 
were AvCoVs most similar to IBV, and were provisionally designated as PhCoV. 

3.2. The two PhCoVs were genetically similar to those of IBVs 

The complete genome sequences of the isolates [0623/17 and [0710/17 were 
obtained by assembly of the fragment sequences. The genome sizes of the two PhCoV 
isolates were both 27,640 bases. Their genome organizations were typical of AvCoVs, 
with the gene order (5' to 3’) of replicase ORF lab, S, E, M, and N (Fig. 2A). Both the 
5‘ and 3’ ends contained short untranslated regions. The replicase ORF lab in both 
isolates occupied 19,892 bp of the genomes (Table 2). This ORF encodes a number of 
putative proteins, including nsp3, nsp5, nsp12, nsp13, and other proteins with 
unknown functions. 

The two PhCoV isolates displayed similar genome organizations to IBV strains, 
and differed only in the number of ORFs downstream of 5a (Fig. 2A). A mutation 
(ATG—GTG) occurred at a corresponding position to the start codon in our two 
PhCoV isolates, resulting in the viruses lacking ORF5b. Another ORF, provisionally 
designated as ORF6, was found between 5a and N in the genomes of the two PhCoVs, 
which overlapped with N. In addition, a 90-bp insertion sequence was found in the 
Sa—ORF6 junction region (Fig. 2B). The genomes of the PhCoV isolates were 
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therefore organized in a novel manner: 5'-Gene 1-S-3a-3b-E-M-4b-5a-6-N-3', which 
differed from those in IBV and other AVCoVs (Fig. 2A). 

In addition, recombination within the genomes of our two PhCoVs was analyzed 
using SimPlot by comparing them with 40 represented IBVs selected on the basis of 
the results of phylogenetic analysis. No obvious recombination events were found in 
the genomes of the PhCoV isolates [0623/17 and I0710/17 (Fig. 2C). 

3.3. The PhCoVs differed from IBV to a similar extent as one lineage of IBV differed 
from another 

The maximum likelihood trees constructed using the SI nucleotide sequences of 
our two PhCoV isolates and 213 reference AvCoV strains are shown in Fig. 3A. The 
two PhCoV isolates were clustered together and were not closely related to IBV and 
other AVCoV strains. Comparatively, the two PhCoV isolates showed a relatively 
closer relationship with IBV GI-19 lineage strains (<92% and 92% nucleotide and 
amino acid identities, respectively). According to the complete genomic sequences, 
the two PhCoV isolates also clustered together and were separate from other IBV 
strains (Fig. 3B). Isolate 10623/17 (10710/17) showed the highest sequence identity 
with the GI-19 lineage (<91%). Trees based on both the S1 and complete genomic 
sequences, thus, showed that the PhCoVs differed from IBV to a similar extent as one 
lineage/genotype of IBV differed from another. 

Nine sequences of S1 gene fragments from PhCoV strains isolated in the UK 
were available in GenBank and used for comparison with our two PhCoV isolates and 
17 IBV strains isolated in China and Europe. All the PhCoVs were divided into three 
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groups (Fig. 3C). Our two PhCoVs were genetically most closely related to GI-19 
lineage viruses, which were believed to originate from China (Liu and Kong, 2004). 
However, four UK PhCoV strains clustered together and showed a closer relationship 
with IBV GI-14 lineage strains isolated in Europe and Africa. The remaining five UK 
PhCoV strains clustered together in a separate group and were most closely related to 
IBV GI-16 lineage viruses, which were isolated in Europe. Multiple alignments using 
these sequences confirmed these results and revealed that PhCoVs in each group 
showed similarities with different IBV lineages (Supplementary Fig. 2), rather than all 
PhCoVs being most closely related to each other. Interestingly, phylogenetic trees 
using the 3' UTR sequences of our two PhCoV isolates, 16 UK PhCoV strains, and 40 
IBV strains (selected based on the results of phylogenetic analysis of the complete 
genomic sequences) showed that most of the PhCoV strains were clustered in a group, 
but four of which clustered differently (Fig. 3D). 
3.5. The PhCoV was antigenically different from the selected IBVs 

The percentage of antigenic relatedness (r) values determined by VN between the 
isolate 10623/17 and other viruses, including isolate 10710/17 and three IBV strains, 
are given in Fig. 4A. The isolate, 10623/17, showed a strong antigenic relationship 
with the isolate, 10710/17, on the basis of r value (67.3). In contrast, the r values of the 
isolate [0623/17 and the IBV strains used in this study were all <18.2%, 1.e. below the 
50% cutoff value considered to indicate greater antigenic dissimilarity. 

An aliquot of each serum was pooled and tested by VN for PhCoV antibodies 
using 10623/17 virus. Neutralization titers ranging from 84.5 (2°*) to 222.8 (27°) were 


22 


obtained from pheasant serum after the disappearance of clinical signs (Fig. 4B), 
indicating previous exposure to the virus. In contrast, the neutralization antibody was 
negative in serum from pheasants that had not shown any clinical signs. 

3.6. The PhCoV was non-pathogenic to SPF chickens but highly pathogenic to 
commercial pheasants 

No clinical signs were observed in chickens inoculated with different PhCoVs at 
different ages (1, 30, and 45 days), comparable with chickens in the negative control 
groups. In contrast, the clinical signs were observed in all the pheasants from 2 to 14 
dpi with isolate [0710/17 at 45 days of age. The clinical signs, including snicking, 
watery eyes, nasal discharge, and rales. No gross lesions were observed in the five 
chickens inoculated with isolate [0710/17 necropsied at 5 dpi, however, gross kidney 
lesions were observed in three of the five pheasants inoculated with isolate I0710/17. 
The kidney parenchyma was pale, swollen, and mottled, and the kidney tubules and 
urethras were distended with uric acid crystals. Three pheasants died at 6 dpi 
following the clinical signs during the experiment. Gross lesions in the dead birds 
were also mainly confined to the kidneys. Mild tracheitis was also observed in the 
dead birds. No clinical signs were observed in the control pheasants. 

The antibody responses are illustrated in Fig. 4C. Both chickens and pheasants 
inoculated with the PhCoV isolates showed negative ELISA antibody results at 4 dpi. 
No antibody responses were observed in any of the chickens inoculated with PhCoV 
isolates [0623/17 and [0710/17 at 8 dpi at different ages in ELISA tests; however, 
50% of pheasants inoculated with I0710/17 showed seroconversion from 8 dpi. The 
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seroconversion rates were higher in pheasants than in chickens at the corresponding 
time points. No antibody was detected in the negative control chickens or pheasants. 
3.7. The PhCoV had high replication capacity in commercial pheasants 

Shedding of the PhCoV isolate in the respiratory and digestive tracts of both 
chickens and pheasants was determined by virus re-isolation using oropharyngeal and 
cloacal swabs from 9-day-old SPF chicken eggs (Han et al., 2016). Prolonged virus 
replication and higher number of birds shedding virus in the respiratory and digestive 
tracts were observed in pheasants, compared to those of chickens inoculated with 
PhCoV 1I0710/17 at different time points (Fig. 5A). Replication of the challenge 
viruses was also determined at 5 dpi in the tracheas, kidneys, lungs, proventriculus, 
and cecal tonsils from chickens and pheasants (Ma et al., 2019) (Fig. 5B). The 
replication capacity of the PhCoV [0623/17 isolate was detected in the tracheas of 
chickens and pheasants with comparable titers; however, the titers were somewhat 
lower in chickens inoculated with I0710/17. Virus titers in the lungs, proventriculus, 
and cecal tonsils were significantly lower in chickens than in pheasants, with no 
viruses detected in some chicken tissues. Notably, no virus was detected in the 
kidneys of chickens challenged with the two PhCoV isolates, in contrast to the high 
titers in the kidneys of pheasants challenged with I0710/17. In line with this result, 
viral antigens were detected by IHC in the kidneys of pheasants challenged with 
10710/17 (Fig. 5C), but not in the kidneys of chickens inoculated with the two PhCoV 
isolates. The viral antigens were also detected by IHC in the tracheas of chickens 
challenged with 10623/17 and pheasants challenged with [0710/17 (Fig. 5C). No virus 


24 


was detected in any tissues or in oropharyngeal and cloacal swabs from birds in the 


control group. 


4. Discussion 

In the current study, we isolated two AVCoVs from commercial pheasants at two 
respective farms in China in 2017 with disease outbreaks occurring within a 2-month 
period. The diseased birds showed early respiratory signs, suggesting that they may 
have been caused by respiratory pathogens, possibly including NDV, AIV subtypes 
H5 and H9, IBV, or PhCoVs. The embryopathic effects and morphologies of the 
isolated viruses and the sequence characteristics of the amplified gene fragments 
suggested that both agents were AvCoVs most similar to IBV. This was confirmed by 
serological ELISA evidence using IBV as the antigen and VN test using the isolate 
10631/17. Since the original description of PhCoV in 1980 (Spackman and Cameron, 
1983), the disease caused by PhCoV has become widespread in the UK (Lister et al., 
1985; Pennycott, 2000; Cavanagh et al., 2002). Some viruses isolated from pheasants 
with respiratory disease, reduced egg production, renal problems, and high mortality 
were believed to be IBV (Spackman &Cameron, 1983; Gough et al., 1996; Lister et 
al., 1985; Pennycott, 2000), due to similarities with IBV in domestic fowl (Cook et 
al., 2012). The mortality associated with the outbreaks at pheasant farms in China 
were comparable to or even higher than those of the outbreak at a large game farm in 
the UK in 1994, which resulted in the deaths of over 1000 breeding pheasants out of a 
total stock of approximately 7000 birds (Gough et al., 1996). The clinical signs and 
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gross lesions in dead birds were also similar between the two outbreaks, with upper 
respiratory signs at the beginning of the outbreaks, and the birds dying rapidly after 
appearance of the signs of illness. Pale, swollen kidneys and ureters distended with 
urates were also found in both cases. Some viruses have been successfully isolated 
and propagated in embryonating domestic fowl eggs (Lister et al., 1985; Gough et al., 
1996), as with the two PhCoVs isolated in this study. 

Phylogenetic analysis showed that our PhCoVs represented a novel lineage in 
genotype I of IBV, according to the method of Valastro et al. (2016). The phylogenetic 
analyses were based on the SI gene sequences of our PhCoVs and 213 reference 
AvCoV strains. The nucleotide and putative amino acid compositions of the SI gene 
also showed approximately 10% differences between the two current PhCoVs isolates 
and other AvCoV strains, compared with 13% of nucleotides and 14% of amino acids 
in the S1 subdomain between IBV lineages, and 30% of nucleotides and 31% 
(occasionally even 40% or more) of amino acids in the S1 subdomain between IBV 
serotypes (Valastro et al., 2016). Our PhCoVs thus differed from IBV to a similar 
extent as one lineage of IBV differed from another. We also found interesting 
branching of the phylogenetic tree based on the complete genomic sequences of our 
two PhCoVs and IBVs of different lineages/genotypes. Pairwise comparisons of the 
complete genomic sequences revealed that our PhCoVs differed by >10%, similar to 
many IBV lineages and serotypes, further supporting the similarity of the genetic 
relationship between our two PhCoVs and IBV strains was similar to that between 
IBV lineages and genotypes. In addition, the robust IB V-specific serum IgG titers in 


26 


chickens inoculated with PhCoV isolates showed that PhCoV infection stimulated a 
humoral response in chickens (Raj and Jones, 1997; Aston et al., 2018). However, 
PhCoVs were antigenically dissimilar to IBVs in the VN tests, suggesting that the 
antigenic relationship between our two PhCoVs and the IBV strains was similar to 
those between different IBV genotypes (Gao et al., 2016). Viruses antigenically 
distinct from the investigated IBVs have also been isolated from pheasants in the UK 
on several occasions (Lister et al., 1985; Cavanagh et al., 2002). 

CoVs typically have a restricted host range and only infect their natural host. As 
an AvCoV, IBV is considered to be a highly contagious virus that only infects 
chickens (Gallus gallus) (Cavanagh, 2007); however, IBV has also been reported in 
other avian species including guinea fowl, partridge, peafowl, and teal, without 
producing any clinical detectable disease (Cavanagh, 2007; Liu et al., 2005; Sun et al., 
2007). Pheasants (P. colchicus) differ from chickens (G. gallus), though both belong 
to the family Phasianidae, order Galliformes. Although the PhCoVs in this study 
were very similar to the IBVs, there were some obvious differences. The IBV RNA 
genome is characterized by the organization 5'-UTR-1a-1b-S-3a-3b-E-M-4b-4c-5a- 
5b-N-UTR-3’, while the genomes of the PhCoV isolates were organized in a novel 5'- 
UTR-1a-1b-S-3a-3b-E-M-4b-5a-6-N-3'UTR order, which differed from those of IBV 
and other AvCoVs. Novel genomic organizations have also been found in other IBVs 
(Liu et al., 2008; Mardani et al., 2008; Xu et al., 2016) and experiments with 
engineered mutants suggested that the native order was not functionally essential 
(Casais et al., 2005, Hodgson et al., 2006). Notably, the current study found 
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differences in the pathogenicity and replication capacity of the PhCoV isolates 
between chickens and pheasants. Intranasal instillation of 10 3-week-old SPF 
chickens with 0.2 ml of infectious allantoic fluid containing 10° EIDs0 of PhCoV 
previously failed to produce any clinical signs of disease (Lister et al., 1985). 
However, all birds had HI, but precipitating antibodies to IBV M41 strain. In this 
study, inoculation of chickens via the ocular and nasal routes with 0.1 ml of 10°° 
EIDs0 of the PhCoV isolates 10623/17 and I0710/17 at 1, 30, and 45 days of age, 
respectively, also failed to produce any clinical signs of disease, while pheasants 
inoculated with I0710/17 showed obvious clinical signs and 30% mortality, similar to 
chickens infected with IBV strains (Liu and Kong, 2004). The clinical signs were 
similar to those in the affected commercial pheasants from which the viruses were 
isolated, and to those in chickens infected with IBV strains (Cavanagh, 2007). 
Furthermore, not all the inoculated chickens had ELISA-detected antibodies, 
compared with seroconversion of all pheasants inoculated with the PhCoV isolate. 
Strikingly, the viruses could be re-isolated from the tracheas of chickens inoculated 
with PhCoV isolates 10623/17 and I0710/17, but were absent or present at very low 
titers in other chicken tissues investigated in this study. In contrast, the PhCoV isolate 
10710/17 could be re-isolated from, and had significantly higher viral titers in all 
pheasant tissues investigated in this study, especially the kidneys. These results were 
in line with the IHC results, which showed that viral antigen-positive epithelial cells 
were intensely labeled in the kidneys of pheasants challenged with I0710/17. These 
results were also consistent with those for viral replication and shedding, which 
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indicated prolonged virus replication and higher number of birds shedding virus in the 
respiratory and digestive tracts of pheasants compared with chickens inoculated with 
PhCoVs. These findings support the different biological characteristics of PhCoVs 
compared with IBVs (Liu et al., 2008; 2011). 

The close relationship—but with obvious differences between PhCoVs and 
IB Vs—led us to speculate that the PhCoVs may have originated from IBVs. This 
hypothesis was confirmed, at least in part, by the analyses of S1 fragments from 
PhCoVs isolated from different outbreaks and selected IBV strains. Instead of being 
clustered together and showing a close genetic relationship, the PhCoVs isolated from 
different outbreaks in different years clustered in different groups, with each group 
showing a relatively close relationship with the IBVs circulating in the regions from 
which the respective PhCoVs were isolated. We speculated that the IBV strains may 
have been introduced into and adapted themselves in pheasants, and subsequently 
evolved independently to different extents resulting in different virulences in 
pheasants. For example, a previous report on the isolation of a PhCoV reported that 
the only notable clinical features were respiratory signs, reduced egg production, and 
shell-quality problems (Spackman and Cameron, 1983). However, the PhCoVs in the 
current and other studies showed a high affinity to the kidney (Gough et al., 1996; 
Lister et al., 1985; Pennycott, 2000) and were highly pathogenic in pheasants. The 
hypothesis that different introductions of the viruses from chickens to pheasants and 
the independent evolution of the resulting PhCoVs may account for the inconsistent 
phenotypes of the PhCoVs. 
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The emergence of CoVs in turkeys was proposed to have resulted from 
recombination events involving IBVs and an as-yet-unidentified CoV, which donated 
an S gene encoding a protein with low amino acid identity to IBV (35%) (Jackwood 
et al., 2010; Brown et al., 2016). This is suspected to have resulted in a host shift from 
chickens to turkeys, as well as switching the tissue tropism of the virus from the upper 
respiratory to the intestinal tract. In this study, we detected genome-wide differences, 
including point mutations, deletions, and insertions, in the ORFs and deduced protein 
sequences of our PhCoVs, compared with IBVs. The accumulation of substitutions in 
multiple genes in the PhCoV genome resulted in a host shift from chickens to 
pheasants, but did not lead to a change in tissue tropism (nephropathogenicity). In 
addition, we found no recombination events in the genomes of our PhCoVs compared 
with the IBV strains selected in this study. In contrast, different phylogenetic 
topologies were observed when the S1 fragments and 3'UTRs of some of the PhCoVs 
isolated in UK were compared, suggesting that possible recombination events might 
have occurred. The different evolutionary patterns between our PhCoVs and those 
isolated in the UK might suggest that PhCoVs isolated in different outbreaks derived 
differently and have evolved independently. 

There is currently no information on the complete genome sequences of PhCoVs 
isolated in other countries. Further molecular studies are therefore required to 
determine the complete genomic sequences of PhCoVs from different outbreaks in 
pheasants, to allow a thorough investigation of the relationships between PhCoVs and 
IBVs. In addition, further research involving the infection of pheasants with IBVs is 
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also required to ascertain its pathogenicity in pheasants and whether could pheasants 


be the potential asymptomatic vectors of IBV for possibly spreading the virus. 
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Figure legends 


Fig. 1. Antibody responses of commercial pheasants, using an ELISA kit for detecting 
IBV antibodies. Blood samples were collected from six male (A and C) and six layer 
hen pheasants (B and D) from each of two affected commercial pheasant flocks in 
Shandong and Anhui provinces, respectively. Blood samples were also collected from 
six additional commercial male (E) and six layer pheasants (F) from other flocks that 
had not shown any clinical signs. Positive and negative serum controls for IBV are 


also shown. 
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Fig. 2. Complete genomic sequence analysis of PhCoV isolate 10623/17 (10710/17) 
with those of other AVCoVs. Genomic organization of PhCoV isolate [0623/17 
(10710/17) compared with AVCoVs TCoV MG1O0, IBV Beaudette, guinea fowl 
coronavirus (GfCoV) Fr2011, and Canada goose coronavirus (CGCoV) 
Cambridge_Bay_2017 (A). The names of the main virion genes are given. Genomes 
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of five AVCoVs are drawn to different scales. Alignment of sequences encoding the 
ORF5 and flanking sequences of PhCoVs 10623/17, 10710/17 and IBV strains 
LDL/091022 AND 1148-A (B). A single mutation at the first site of the start codon 
corresponding to the IBV 5b (ATG—ATT) resulted in absence of 5b in the PhCoV 
genome. A 90-nt insertion was found in the 5a-ORF6 junction region in the PhCoVs. 
Analysis was also performed using SimPlot software version 3.5.1 to identify 
potential recombination breakpoints in the genomes of PhCoV isolates 10623/17 and 
10710/17 (C). We used a 1000-bp window with a 50-bp step. Isolate I0710/17 was 
used as the query strain, and 39 selected IBV strains and one duck CoV were used as 


reference strains. The names of the main virion genes are given at the top of the 


SimPlot. 
Fig. 2A. 
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Fig. 3. Phylogenetic analysis. Two hundreds and thirteen complete S1 nucleotide 
sequences of AVCoV reference strains were phylogenetically compared with those of 
PhCoV isolates [0623/17 and [0710/17 (A). The strains LDL/091022 and LJL/140734 
(star) were used in the following VN tests. The complete genomic sequences of our 
two isolates 10623/17 and I0710/17 were also phylogenetically compared with those 
of 85 AvCoV reference strains, including 84 IBV strains and one duck CoV strain (B). 
The 40 sequences used for SimPlot analysis are shown in bold. Phylogenetic trees 
were constructed using the nucleotide sequences of the S1 gene fragments of PhCoV 
isolates 10623/17 and I0710/17 (black circle), nine PhCoV strains isolated in the UK 
(triangle), and 11 selected IBV strains (C). Phylogenetic trees were also conducted 
using the sequences of 3’ UTR fragments from PhCoV isolates [0623/17 and [0710/17 
(black circle), 16 PhCoVs isolated in the UK (triangle), 43 selected IBV strains, and 
one duck CoV (D). The maximum likelihood method was used to construct the trees 
with 1000 bootstrap replicates with MEGA4.0 software. GenBank accession numbers 
are indicated after the names of the viruses PhCoV isolates 10623/17 and 10710/17 


(black circle) 
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Fig. 4. Neutralization and ELISA tests. The neutralization tests were conducted 
between PhCoV isolates 10623/17 and I0710/17 and between isolate 10623/17 and 
selected IBV strains (A). The calculated antigenic relatedness (r) value of PhCoV 
isolate [0623/17 against I0710/17 and selected IBV strains were shown. The 
neutralization titers of serum from the commercial pheasants tested by VN for PhCoV 
antibodies using 0623/17 virus were tested (B). Sera ELISA antibody responses of 
chickens inoculated with PhCoV isolates [0623/17 and [0710/17 at different ages and 


pheasants inoculated with PhCoV isolate 10710/17 at 45 days of ages (C). 
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Fig. 4B. 
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Fig. 5. Infection of chickens with PhCoV isolates 10623/17 and I0710/17 and 
pheasants with I0710/17. Virus recovery from oropharyngeal (I) and cloacal (II) 
swabs from chickens inoculated with PhCoV isolates [0623/17 and I0710/17, and 
pheasants with PhCoV isolate I0710/17 (A). Virus recovery was performed by 
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inoculating 9-day-old embryonated, specific-pathogen-free eggs through the allantoic 
route with supernatant from the swabs. Replication of our isolates 10623/17 and 
[0710/17 in the trachea, lung, proventriculus, cecal tonsil, and kidney of chickens at 5 
dpi, and replication of I0710/17 in the trachea, lung, proventriculus, cecal tonsil, and 
kidney of pheasants at 5 dpi (B). Viral titration was performed by inoculating 9-day- 
old embryonated, specific-pathogen-free eggs through the allantoic route. Data are 
expressed as the mean + standard deviation. Virus titers were analyzed by Student’s f- 
test using GraphPad Prism for Windows version 5. Differences were considered 
significant if the p value was <0.05 (*p<0.05, **p<0.01, ***p<0.001). Detection of 
the IBV antigen by IHC analysis in the kidney (J) and trachea (II) of pheasants at 5 
dpi following PhCoV [0710/17 infection and in the trachea (III) of chickens at 5 dpi 


following PhCoV [0623/17 infection (C). 
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Fig. 5C. 


50 


Table 1. 


Information of experimental design. 


Group —_Chicken/Pheasant Number® Age? Challenge virus 

1 Chicken 15 1 10623/17 

2 Chicken 15 1 10710/17 

3 Chicken 15 ] Sterilized allantoic fluid 
4 Chicken 15 30 10623/17 

5 Chicken 15 30 10710/17 

6 Chicken 15 45 Sterilized allantoic fluid 
dq Chicken 15 45 10623/17 

8 Chicken 15 30 10710/17 

9 Chicken 15 30 Sterilized allantoic fluid 
10 Pheasant 15 45 10710/17 

1] Pheasant hs 45 Sterilized allantoic fluid 


A total of 15 birds were included in each group. Five birds from each group were 
euthanized at 5 dpi, and trachea, lung, proventriculus, kidney, and cecal tonsil tissues 
were collected and used for virus titration. 


> The age when the birds were challenged with the virus. 
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Table 2. 
Genome organization and predicted viral proteins encoded by pheasant coronaviruses 


10623/17 (10710/17). 


Gene Genome Size Size (amino Identity to|IBV  Most-related 
position (nucleotide) acid) (%) lineage/strain 
5'UTR 1-529 529 - 97.9 GI-19 
ORF NSP2 530-2551 2022 674 90.0 GI-19 
1a/1ab NSP3 = =2552-7318 4767 1588 90.6 GI-22 
NSP4 7319-8866 1548 516 92.7 Gl-19 
NSP5 8867-9787 921 307 co GI-22 
NSP6 9788-10669 882 284 94.7 Gl-22 
NSP7 10670-10918 249 83 98.8 1114/14 
NSP8 10919-11548 630 210 93.2 Gl-22 
NSP9 11549-11881 333 111 92.8 DK/CH/HN/ZZ2004 
NSP1 11882-12316 435 145 94.5 DK/CH/HN/ZZ2004 
0 
NSP1 12317-12385 69 23 7 7 
1 
NSP1 12317-15135 2819 931 94.9 GI-19 
2 
NSP1 15136-16935 1800 600 93.0 GI-19 
3 
NSP1 16936-18498 1563 521 93.9 Gl-22 
4 
NSP1 18499-19512 1014 338 93.1 GI-19 
5 
NSP1 19612-20446 834 278 95.2 GI-19 
6 
S S1 20372-21991 1620 540 92.4 GI-19 
S2 21992-23869 1878 626 98.4 GI-19 
ORF 3 3a 23869-24042 174 58 99.4 GI-19 
3b 24042-24236 195 65 100 GI-9 
3c(E) 24217-24546 330 110 95.7 GI-19 
M 24518-25195 678 226 95:4 GI-29 
ORF 4b 25196-25480 285 95 91.2 GI-19 
ORF 5a 25559-25756 198 66 92.4 GI-19 
ORF 6 25849-26091 243 81 87.0 Gl-19 
N 26034-27260 1227 409 92.5 GI-19 
3' UTR 27261-27640 380 - 98.4 GI-19 
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